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order including all relevant dimension 4 current corrections. 
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Table 1: Simulations details on three "coarse" and three "fine" MILC ensembles. 



1. Introduction 

Leptonic decays of charged B's such as — )• t+Vt, are important processes for Cabbibo- 
Kobayashi-Maskawa (CKM) and Unitarity Triangle (UT) physics. There is currently some tension 
between Ek, sin(2j8), \Vub\ and ^(B+ — )• T+v^) and the B meson decay constant /g plays an impor- 
tant role in such analyses. For instance global fit results for fs from precision electroweak data are 
being compared with "direct Standard Model predictions" of the decay constant from lattice QCD 
[|l|, ^ ^. Reducing errors in the lattice determinations of fs is, hence, a worthwhile high priority 
goal. 

For the meson there are no tree-level leptonic decays in the Standard Model (SM). Never- 
theless the decay constant /b, is a useful parameter in many decay and mixing rates both within 
and beyond the SM. Furthermore, lattice QCD determinations of /g, can be achieved with smaller 
errors than for fs since no light valence quarks are involved here. This combined with the fact that 
the ratio is also known more accurately (due to cancellation of statistical and many systematic 
errors) than either of the decay conatants on their own, opens the possibility of getting at precision 
fs via precision values for fs^ and 

The HPQCD collaboration has initiated new calculations of the B and B, meson decay con- 
stants based on NRQCD bottom and HISQ light and strange valence quarks and employing MILC 
Nf = 2 + \ configurations. This improves on previous determinations of fg, /b, and by HPQCD 
that used AsqTad light and strange valence quarks [Q]. Simulation details are summarized in Table 
1. 

2. Tuning of Quark Masses 

We use the static potential quantity ri = 0.3133(23)fm |^ to set absolute scales and ri/a from 
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Figure 1: Tuning of the b-qnark mass using the spin averaged T mass. 

MILC 1^ for relative scales between different MILC ensembles. To fix the bare ft-quark mass in 
lattice units aMi, we use the spin averaged T mass. One calculates, 

^hb = \ [^MtinCSi) +MkinCSo)] , (2.1) 

with 

Mm = ' , AEp = E{p)-E{0), (2.2) 

and compares with the experimental value (adjusted for the absence of electromagnetic, annihila- 
tion and sea charm quark effects in our simulations) of 9.450(4)GeV Results from this tuning 
are shown in Fig. 1. Errors in the data points include statistical and ri/a errors. One sees that 
these are much smaller than the 0.7% error in the absolute physical value of ri. To achieve small 
statistical errors in Mm it was crucial to employ random wall sources for the NRQCD Z^-quark 
propagators. With point sources errors would have been about 4~5 times larger The 5-quark mass 
was tuned to the (fictitious) rj^ mass of 0.6858(40)GeV Having fixed the bottom and strange 
quark masses on each ensemble one can check how well Mb^ —Mj^-j^/l is reproduced in the contin- 
uum limit. The leading dependence on the heavy quark mass cancels in this difference, so one is 
testing how well the lattice actions are simulating QCD boundstate dynamics. Results for this mass 
difference are shown in Fig. 2. Within the ri scale error and additional ~10MeV uncertainty from 
relativistic corrections to M^^ one sees agreement with experiment in the continuum limit. 

3. The Currents and Matching 

In the Bq meson rest frame (q = light or strange) the decay constant is defined in terms of the 
temporal component of the bq heavy-light axial vector current Aq as, 

{0\Ao\B,)qcd=MbJb,- (3.1) 



> 
O 

3 



9.5 



■T3 

<u 

> 
< 

e 

00 



9.4 



n a 



3 



Studies ofB and Leptonic Decays with NRQCD and HISQ Quarks 



Junko Shigemitsu 



0.68 r 



> 0.67 - 

o 



l;g 0.66 



I 0.65 



J Coarse 



Fine 
experiment 
0.7% Scale Error 



# 



pa 



0.64"! 



0.05 



0.1 



0.15 0.2 
(a/r.)**2 



0.25 



Figure 2: The mass difference Mb, —Myj^/l versus the square of the lattice spacing for the first five ensem- 
bles of Table 1 . One sees negligible sea quark mass dependence but a noticeable lattice spacing dependence. 



Simulations are carried out with effective lattice theory currents, 

J^^\x) = q{x)ToQ{x), (3.2) 

= ^^Wror-veW, (3.3) 

4"' (^) = :^^(^) r- ^ ro To Q{x) , (3.4) 
ZMh 

and matching through order a,, aa^, Ois^^^ gives, 

(Ao)ecD = (1 + asPo) (/f) + (1 + «.Pi) + «^P2 (4'^'™'). (3-5) 



i(i);sub _ 



«.Cio/(°^. (3.6) 



Po> Pi» P2 and (^10 are the one-loop matching coefficients which have recently been calculated for 
NRQCD/HISQ currents. 

4. Preliminary Results and Error Estimates 

Figs. 3 & 4 show our preliminary chiral/continuum extrapolations of /b, ^/Mb,, /bV^b and the 
ratio /b, \/MbJ /bV^b based on the first five ensembles of Table 1. Work on the sixth ensemble 
FO, a more chiral fine ensemble, is in progress. Extrapolations are carried out using continuum 
partially quenched ChPT for heavy-light decay constants augmented by lattice spacing dependent 
terms. Our preliminary numbers at the physical point are, 

/B = 0.191(9)GeV, /B,=0.226(10)GeV, = 1.184(19). (4.1) 
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Figure 3: Continuum/chiral extrapolations of fs, \JMb^ (upper curves) and Jb^JMb (lower curves). The 
black circles are results at the physical point from reference |Q| using AsqTad light and strange valence 
quarks. 
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Figure 4: Continuum/chiral extrapolation of the ratio /b, ^Mg, / /b \/Mb ■ The black circle is the result at 
the physical point from reference M] using AsqTad light and strange valence quarks. 
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Table 2: Preliminary error budget 



Table 2 shows a preliminary error budget. 

The new NRQCD/HISQ numbers are consistent with HPQCD's NRQCD/AsqTad results of 
/b = 0.190(13)GeV, fe, = 0.231(15)GeV, /bJ/b = 1.226(26) !§. The reduction in errors in the 
new calculations comes mainly from improvement in discretization errors, smaller scale un- 
certainties and better fitting and extrapolation strategies. The striking decrease in lattice spacing 
dependence as one goes from AsqTad to HISQ strange quarks is demonstrated in Fig.5, where 
we compare fs^ \/Mb, results on the same ensembles for the two different valence quark actions. 
Although this improvement in discretization errors is very welcome, in both the present and pre- 
vious calculations the total error is dominated by the higher order operator matching uncertainty. 
HPQCD is investigating nonperturbative matching strategies for NRQCD/HISQ currents which 
could reduce this error in the future 

In a completely different thrust, we are pursuing an alternate approach to B physics that uses 
the relativistic HISQ action for heavy quarks with masses niH > mcharm on very fine lattices [^, p^ ]. 
One can then extrapolate up to the physical ^-quark staying always within a relativistic frame 
work. A recent result using this method gives /b^ = 0.225 (4)GeV [ pl| ] in excellent agreement with 
eq.(4.1), however with much reduced errors. Using relativistic heavy quarks enables us to work 
with absolutely normalized currents (based on Ward identities). The main source of the larger 
errors with NRQCD Z7-quarks, namely operator matching uncertainties, is thus removed. In order 
to carry out heavy HISQ calculations that can be extrapolated to the physical b-quaik, very fine and 
hence large lattices are required. Repeating the /b^ determination for /g would be quite expensive. 
So for the next couple of years we believe the best strategy for precision B physics will be to 
work with HISQ heavy quarks for Bs physics and combine these results with ratios, such as /bJ/b 
from NRQCD ft-quark calculations. The calculations presented here are part of this comprehensive 
approach to precision B physics. 
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Figure 5: Comparison of lattice spacing dependence between NRQCD/HISQ results for /bs^/Mb, from 
the current calculations (open circles and triangles) and previous results from reference based on 
NRQCD/AsqTad quarks. 



References 

[1] E. Lunghi, A. Soni, Phys. Lett. B697, 323-328 (2011). [arXiv: 1010.6069 [hep-ph]]. 

[2] E. Lunghi, A. Soni, Phys. Rev. Lett. 104, 251802 (2010). [arXiv: 09 12.0002 [hep-ph]]. 

[3] J. Laiho, E. Lunghi, R. Van De Water, PoS FPCP2010, 040 (2010). [arXiv: 1 102.3917 [hep-ph]]. 

[4] E. Gamiz et al. [ HPQCD Collaboration ], Phys. Rev. D80, 014503 (2009). [arXiv:0902.1815 
[hep-lat]]. 

[5] C. T. H. Davies et al. [ HPQCD Collaboration ], Phys. Rev. D81, 034506 (2010). [arXiv:091 0.1229 
[hep-lat]]. 

[6] A. Bazavov et al. [ MILC Collaboration ], Rev. Mod. Phys. 82, 1349-1417 (2010). [arXiv:0903.3598 
[hep-lat]]. 

[7] E. B. Gregory et al. [ HPQCD Collaboration ], Phys. Rev. D83, 014506 (201 1). [arXiv: 1010.3848 
[hep-lat]]. 

[8] J. Koponen et al. [ HPQCD Collaboration ], PoS LATTICE2010, 23 1 (2010), [arXiv: 101 1 . 1208 
[hep-lat]] and K. Hornbostel et al. [ HPQCD Collaboration ] work in progress. 

[9] C. McNeile et al. [ HPQCD Collaboration ], PoS LAT2009, 1 16 (2009). [arXiv:0910.2921 [hep-lat]]. 

[10] E. Follana et al. [ HPQCD Collaboration ], PoS LAT2010, 305 (2010). 

[11] C. McNeile et al. [ HPQCD Collaboration ], arXiv: 1 1 10.4510 [hep-lat]. 



7 



